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Echoes from the companion star in Sco X-l 
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ABSTRACT 

We present simultaneous X-ray (RXTE) and optical (ULTRACAM) narrow band (Bowen 
blend/Hell and nearby continuum) observations of Sco X-l at 2-10 Hz time resolution. We 
find that the Bowen/Hell emission lags the X-ray light-curves with a light travel time of ~ 
11 — 16s which is consistent with reprocessing in the companion star. The echo from the 
donor is detected at orbital phase ~ 0.5 when Sco X-l is at the top of the Flaring Branch. 
Evidence of echoes is also seen at the bottom of the Flaring Branch but with time-lags of 5-10 
s which are consistent with reprocessing in an accretion disc with a radial temperature profile. 
We discuss the implication of our results for the orbital parameters of Sco X-l. 
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1 INTRODUCTION 

Low mass X-ray binaries (LMXBs) are interacting binaries con- 
taining a low mass donor star transferring matter onto a neutron 
star (NS) or black hole. Mass accretion takes place through an ac- 
cretion disc, and with temperatures approaching ~ 10 8 K, such 
systems are strong X-ray sources. The mass transfer rate sup- 
plied by the donor star, M2, is driven by the binary /donor evo- 
lution and for Mi > Merit ~ lO _9 M0j/r _1 results in per- 
sistently bright X-ray sources. In these binaries optical emission 
is dominated by reprocessing of the powerful X-ray luminosity 
in the gas around the compact object which usually swamp the 
spectroscopic features of the weak companions stars. Only in ex- 
ceptional cases, such as long-period LMXBs with much brighter 
evol ved donors (e. g. Cyg X-2) is a com plete orbital solution possi- 
ble dCasares. Charles & Kuulkers 19981 : hereafter CCK98). In this 
scenario, dynamical studies have classically been restricted to the 
analysis of X-ray transients during quiescence, where the intrinsic 
lumi nosity of the donor do minates the light spectrum of the binary 
(e.g. lCharles&Coe 2004b . 
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1.1 Fluorescence Emission from Donor Stars 

Fortunately, this situation has recently changed thanks to the 
discovery of the narrow emission components arising from 
the donor star in the prototypical persi stent LMXB Sco X-l 
IfSteeghs & Casares 2002; hereafter SC02I High resolution spec- 
troscopy revealed many narrow high-excitation emission lines, the 
most prominent associated with Hell A4686 and NIII AA4634-41 / 
CIII AA4647-50 at the core of the broad Bowen blend. In particular, 
the NIII lines are powered by fluorescence resonance through cas- 
cade recombination which initially requires EUV seed photons of 
Hell Ly«. These very narrow (FWHM < 50km s _1 ) components 
move in antiphase with respect to the wings of Hell A4686, which 
approximately trace the motion of the compact star. Both properties 
(narrowness and phase offset) imply that these components origi- 
nate in the irradiated face of the donor star. This work represented 
the first detection of the companion star in Sco X-l and opened 
a new window for extracting dynamical information and thereby 
deriving mass functions in a population of ~ 20 LMXBs with es- 
tablished optical counterparts. We now know that this property is 
not peculiar to Sco X-l but is a feature of persistent LMXBs, as 
demonstrated by the following examples: 

(i) Doppler tomography of NIII A4640 in the eclipsing ADC 
(accretion disc corona) pulsar 2A 1822-371 reveals a com- 
pact spot at the position of the comp anion star and wit h 
velocity K em = 300 ± 8 km s" 1 dCasares et ai] l2003h . 
Moreover, the application of the K-correction to this K em 
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velocity, combined with other system parameters, strongly 
points to the presence of a neutron star with mass 1.6 — 
2.3Mr in this system jMunoz-Darias. Casares. & Martinez-Paid 
12005; hereafter MCMOSri . 

(ii) Radial velocities of narrow Bowen lines in the black hole 
candidate GX339-4, detected during the 2002 outburst, led to 
a mass function in excess of 5.8Mq, and hence provided the 
first dynamical proof that the compact object is indeed a black 
hole dHvnes et al. 20031) . 

(iii) Sharp NIII A4640 Bowen emission has bee n detected 
in 4U 1636-536 4U 1735-444 dCasares et al] 120061) . Aql X-l 
dCornelisse et al.l 20071) and th e transient millisecond pulsar XTE 
J1814-338( Casares et al. 20041) . which lead to donor velocity semi- 
amplitudes in the range 200-300 km s _1 . 



Table 1. Sco X- 1 optical observing Log 



1.2 Echo-Tomography 

Echo-tomography is an indirect imaging technique which uses 
time delays between X-ray and UV/optical light-curves as a 
function of orbital phase in o rder to map the reprocessing sites in 
a binary dO'Brien et al. 200*21 hereafter OB02). The optical light 
curve can be simulated by the convolution of the (source) X-ray 
light curve with a transfer function which encodes information 
about the geometry and visibility of the reprocessing regions. 
The transfer function (TF) quantifies the binary response to the 
irradiated flux as a function of the lag time and it is expected 
to have two main components: the accretion disc and the donor 
star. The latter is strongly dependent on the inclination angle, 
binary separation and mass ratio and, therefore, can be used to 
set tight constraints on these fundamental parameters. Successful 
echo-tomography experiments have been performed on several 
X-ray active LMXBs using X-ray and broad-band UV/optical 
light-curves. The results indicate that the reprocessed flux is 
dominated by the large contribution of the accretion disc (e.g. 
Hvnes et al 1998l :hereafter H98, OB02, lO'Brien et al. 20041 and 
Hvnes 2005) which di lutes the reproc e ssed s ignal arising from the 
compa nion. Only iMiddleditch et al. (198lj) and iDavidsen et all 
( 1975) have reported the detection of reprocessed pulsations in 
the X-ray binaries 4U 1626-67 and Her X-l respectively. They 
only used optical data but the Doppler phase modulation of the 
coherent pulses leads to reproce ssing of the X-ray pulses on the 
companion star. More recently, Hvn es et al.l j2006h have reported 
some evidence for a reprocessed signal from the companion 
through the detection of orbital phase dependent echoes in EXO 
0748-676. 

Exploiting emission-line reprocessing rather than broad-band 
photometry has two potential benefits: a) it amplifies the response 
of the donor's contribution by suppressing most of the background 
continuum light (which is associated with the disc); b) since the 
lines are formed in an optically thin medium, the reprocessing 
time is almost instantaneous and hence the response is sharper 
(i.e. only sme ared by geometry) and the transfer function easier to 
compute (see lMunoz-Darias et al] |2005l ; hereafter MD05). Taking 
advantage of this, we decided to undertake an echo-tomography 
campaign on the brightest LMXB of all, Sco X-l, in order to 
search for reprocessed signatures of the donor using simultaneous 
X-ray and Bowen/Hell emission line light-curves. 



Date(2004) 


Texp(s) 
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$ orbital 
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17 May 


0.1 
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0.07-0.35 


60 


18 May 


0.25-0.5 


1-3 


0.34-0.73 


50-100 


19 May 


0.3 


1-2 


0.73-0.95 


90 



1.3 Background to Sco X-l / V818 Sco 

Sco X-l is the prototype LMXB and also the brightest persis- 
tent X-ray source in the sk y and has been the t arget of detailed 
studies since its discovery dGiacconi et al. 19621). An orbital pe- 
riod of 18.9 hr was reported by iGottlieb et alj|( 1975)1 from photo- 
graphic photometry of the V~13 optical counterpart, V818 Sco, 
using plates spanning over 85 years. The opti cal B-band light 
curve shows a 0.13 mag amplitude modulation dAugusteiin et al.l 
1 19921) . which is interpreted as due to the changing visibility of 
the inner-hemisphere of the X-ray heated companion. Radial ve- 
locity curves of HeIIA4686 and HI confirmed the orbital period 
and indicated that the inferior conjunction of the emission line 
regions are close to the photometric minimu m and, hence, the y 
must originate near t he compact object (e.g. LaSa la et al. 1985). 
iBradshaw et alj(1999)l measured the trigonometric parallax of Sco 
X-l using VLBA radio observations and deduced a distance of 
2.8 ± 0.3 kpc. Further observations revealed the presence of twin 
radio lobes, whic h are inclined at an angle of 44° ± 6° relative to 
the line of sight dFomalont et al. 200 ll ; hereafter F01). Furthermore, 
Sco X-l is classed as a Z — source due to the Z-shape patterns that 
the spectral variations trace in the X-ray hard/ soft colour-colour 
diagr am. The source also exhibits kHz QPOs dvan der Klis et al.l 
1996) but no X -ray bursts have been observed so far. 
Possible correlated optical and X-ray variability has been stud- 
ied by using broad band observations in both spec tral ranges. 
In particular. Illovaiskv et al] dl980; he reafter 180) and IPetro et all 
( 1981; hereafter P8ll) present evidence for correlated variability 
during flaring episodes. They measured time-lags consistent with 
reprocessing in the accretion disc. Here we present the first detec- 
tion of delayed echoes arising from the companion star in Sco X-l 
together with other echoes more consistent with reprocessing in the 
accretion disc. The echoes are found during different X-ray states 
within the Z-shape pattern. We discuss the implications of our re- 
sults for reprocessing theory and the system parameters of Sco X-l. 



2 OBSERVATIONS AND DATA REDUCTION 

Simultaneous X-ray and optical data of Sco X- 1 were obtained on 
the nights of 17-19 May 2004. The full 18.9 hr orbital period was 
covered in 12 snapshots, yielding 16.1 ks of X-ray data simultane- 
ous with optical photometry in 3 different bands. 



2.1 X-ray data 

The X-ray data were obtained with the Proportional Counter Array 
(PCA) onboard th e Rossi X-ray Timing Explorer (RXTE) satellite 
djahoda et aHll996h . Only 2 PCA detectors (2 and 5) were used 
and the pointing offset was set to 0.71° due to the brightness of Sco 
X-l. The data were analysed using the FTOOLS software and the 
times corrected to the solar bary centre. The STANDARD- 1 mode, 
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Figure 1. Simultaneous RXTE and ULTRACAM (SDSS u', Bowen+Hell 
and Continuum) observations of Sco X-1/V818 Sco obtained on the night 
of 18 May 2004. 

with a time resolution of 0.125s, was used for the variability anal- 
ysis. Additionally, the PCA energy bands 6(2-4 keV), c(4-9 keV) 
and d(9-20 keV) of the STANDARD-2 mode data, with a time res- 
olution of 16s, were used for the spectral analysis. Sco X-l was ob- 
served during 15 RXTE windows of 16-31 minutes length spread 
over the three nights and yielding 20 ks of data. About 80% of these 
observations were also covered with simultaneous optical data. 

2.2 Optical data 

The optical data were obtained with ULTRACAM at the Cassegrain 
focus of the 4.2m William Herschel Telescope (WHT) at La 
Palma. ULTRACAM is a triple-beam CCD camera which uses two 
dichroics to split the light into 3 spectral ranges: Blue (< A3900), 
Green(AA3900-5400) and Red (> A5400). It uses frame trans- 
fer 1024x1024 E2V CCDs which are continuously read out, and 
are capable of time resolution do wn to 2 millisecond s by reading 
only small selected windows (see lDhillon et alj|2007l for details). 
ULTRACAM is equipped with a standard set of ugriz SDSS fil- 
ters. However, since we want to amplify the reprocessed signal 
from the companion, we decided to use two narrow (FWHM =100 
A) interference filters in the Green and Red channels, centered 
at A e ff=4660A and A c ff=6000A. These filters block out most of 
the continuum light and allow us to integrate two selected spec- 
tral regions: the Bowen/Hell blend and a featureless continuum, 
from which continuum-subtracted light-curves of the high excita- 
tion lines can be derived. A standard SDSS u'(3543 A) filter was 
also mounted in the blue channel. The images were reduced using 
the ULTRACAM pipeline software with bias subtraction and flat- 
fielding. Star counts were extracted by adjusting the aperture radii 
to the seeing, and light-curves were obtained relative to a compar- 
ison star which is 96 arcsecs NW of Sco X-l. Light-curves of the 
Bowen/Hell lines (hereafter B+Hell CS) were computed by sub- 
tracting the Red (Continuum) from the Green(hereafter B+Hell) 
channels (see section 3.1.2). The seeing was 1-1.5 arcsecs most of 
the time, except for the first two RXTE windows of 1 8 May when 
it degraded to over 3 arcsecs. Optical observations during the first 3 
RXTE visits on 19 May were ruined by clouds. The exposure time 
was initially set to 0.1s but was increased to 0.25s depending on 
weather conditions. Integrations of 0.5s had to be used for the sec- 
ond window on 18 May, when the seeing was worse. An observing 
log is presented in TableQ] 
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Figure 2. Zoom of W3 showing a 1000s stream of X-ray (above) and 
B+Hell non continuum subtracted (below) data. Both low and high fre- 
quency variability appear clearly correlated. 

3 LIGHT-CURVES 

In order to look for correlated variability we performed a compre- 
hensive cross-correlation analysis by using the entire data set in 
blocks from 2 min to 30 min (i.e. one RXTE window). Significant 
cross-correlation peaks were found only on the night of 18 May, 
when the X-ray flux rose up to 1.3 x 10 4 counts s _1 , which is a 
factor ~ 2 higher than the mean flux displayed by the source dur- 
ing the whole campaign. During this night the system also showed 
long episodes of flaring activity (see figQJ- Unfortunately, during 
the first two RXTE visits the weather conditions were very poor 
and no useful optical data could be obtained. In this work we are 
going to focus on both the 3rd and 5th RXTE visits of May 18th 
where we detect clear episodes of correlated variability between the 
X-ray and optical data. We also note that some correlated variabil- 
ity appears in the 4th RXTE visit but in this case the correlation 
peaks are noisier and many uncorrelated features are also present 
in the data. 

3.1 Echo detection of the companion at superior conjunction 

In fig. [2] we present a zoom of the large flare seen during the 3rd 
RXTE visit of 18 May (W3) together with the simultaneous B+Hell 
light curve. During W3, where a big flare occurs, the X-ray flux 
reached 1.3 x 10 4 counts s _1 , the maximum count rate in the 
whole campaign. We note significant X-ray variability with clear 
correlated counterparts in the optical data. According to the or- 
bital ephemeris reported by SC02, W3 is centered at orbital phase 
0.52 i.e. very close to the superior conjunction of the donor star, 
when the irradiated face presents its largest visibility with respect 
to us. Assuming a light path from the X-ray source to the center of 
the companion, the expected time delay (r) for reprocessing in the 
donor is: 

t — — (1 — sin i cos 2n(j> or b) (1) 

where i is the inclination angle, a the orbital separation, 4>orb the 
orbital phase and c the velocity of light. Therefore, (j>orb ~ 0.5 
is the phase when reprocessing in the companion star is expected 
to be at maximum. For instance, if we use standard values for 
a(~ 10 lt-sec) and i(~ 44°), we obtain r ~ 17s, which should 
be considered as an upper limit to r since the reprocessing take 
place on the heated hemisphere of the donor and not in its center 
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Figure 3. ICF analysis for W3. The thick solid line presents the ICF ob- 
tained between the X-ray and B+Hell emission. We also plot the ICF of 
the low (dashed line) and high frequency (thin solid line) variability after 
filtering the light curve with a step function at 0.002Hz. A clear peak corre- 
sponding to a delay of 10-1 Is (dotted line) is obtained in the latter case. 

of mass. 



3.1.1 Cross-correlation analysis 

As a first step we have computed cross correlation functions be- 
tween the X-ray and the B+HeII light-curves by using the en- 
tire data set obtained during W3. We have used a modified 
version of the ICF (Interpolation Co rrelation Function) method 
which is explained in great detail by Ga skell & Peterson! If 1 987ll . 
We also tested the DCF (Discrete Correlation Function) method 
jEdelson & Krolikll 19881) but found no differences in their results. 
Fig.[3]shows the ICF obtained for W3 as the thick solid line, where 
a clear peak centered at 8-30 s is obtained. Since the light-curves 
in W3 show evidence for both high and low frequency variability 
we decided to filter, prior to computing the ICFs, in order to decon- 
volve the two components.To do this, we convolved the data with 
a sine function to isolate the low frequency component. The light- 
curve corresponding to the high frequency variability was obtained 
by subtracting the low frequency component light-curve from the 
original data. The ICF computed for frequencies lower than 0.002 
Hz, which roughly traces the profile of the big flare, is plotted in 
fig. [3] as the dashed line. This shows a strong and wide peak cen- 
tered at ~ 40 — 80s. This time-lag is clearly longer than the re- 
processing time scale (RTS) which, for Sco X-l, is expected to 
be in the range — 20s. Hence, it cannot be associated with the 
light travel time alone and its origin will be discussed in section 
5. The high frequency (> 0.002 Hz) analysis, on the other hand, 
yields a flat ICF (fig.[3]thin solid line) with a narrow peak centered 
at a lower delay. We have fitted this peak to both parabolic and 
gaussian functions finding a time-lag between 10 — lis. This delay 
is consistent with reprocessing on the companion star (see section 
5). Based on this result we decided to focus our analysis on short 
(~ 2 — 5 min) data intervals dominated by high frequency vari- 
ability. We obtain highly significant correlation peaks by choosing 
different data blocks during W3 especially during the second part 
of this window where strong correlated variability is present(see fig 
[2). In many cases we find ICF time-lags of about lis between the 
X-ray and B+HeII light-curves. For detailed analysis we finally se- 
lected a 3 min block where the correlation peaks are particularly 
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Figure 4. The 3 min interval, selected from W3, for the ICF analysis. 

significant. This interval is presented in fig.E] The resulting ICF 
function is shown in figure[5]as the thin line. A Gaussian fit to the 
peak yields a time lag of 11.8 ± 0.1s, where the error represents the 
formal error of the fit. In order to establish the confidence levels 
we have performed a Monte Carlo simulation where synthetic ICF 
were computed from a large population (10 J ) of white noise time 
series (random, zero-mean normal distributions of data) with the 
same sampling as the optical data. The correlation of each random 
function with respect to the original X-ray data was computed and 
then a statistic was calculated. Thus a 3a level corresponds to the 
upper envelope which encompasses 99.7% of the computed ICFs. 
For the selected interval we find a 3<r confidence at a correlation 
level of 0. 125 which is much lower than our ICF peak value. 
The same analysis was also performed with the light-curves ob- 
tained using the SDSS u and Continuum filters. The correspond- 
ing ICFs also show clear peaks centered at positive delays with 
correlation levels very close to 0.5. However, in these two cases the 
correlation peaks are wider and non-symmetric, and hence the error 
obtained from the Gaussian fits are larger. In figure|5]we also show 
the ICF for the Continuum data which lags the X-ray emission by 
9.0±0.7s. In the same figure we have over-plotted as the dotted line 
the ICF obtained between the Continuum and the B+HeII emission 
where we measure a delay of 2.2 ± 0.3s. This time-lag is consistent 
with the above results. On the other hand, we find that the u band 
emission lags the X-ray data by 9.5 ± 0.7s, which is also ~ 1 — 3s 
shorter than the delay obtained with the B+HeII filter. 



3.1.2 ICF analysis for continuum subtracted B+HeII 
light-curves 

SC02 pointed out that at least 10% of the flux contained in the 
Bowen/Hell emission lines arise from the companion star. Thus, 
in an attempt to amplify the signal from the donor we decided to 
subtract the Red (Continuum) from the Green(B+HeII) channels. 
However, we have seen above that the Continuum is strongly 
correlated with the X-ray emission, and hence, the subtraction 
process will reduce the correlation level of the signal and will add 
noise. Furthermore, the continuum level will likely be different at 
the central wavelength of the red continuum filter and the Bowen 
spectral region. The ratio of continuum levels probably depends 
on the X-ray state of the source as well as perhaps the orbital 
phase since the companion contribution changes around the orbit. 
Therefore, we performed a test to determine the optimal amount 
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Figure 5. Continuum versus Bowen+Hell ICFs for the selected 3min block 
of W3. The solid lines show the correlations between the X-ray and both 
the B+Hell and Continuum emission. As a test we also plot (in dotted line) 
the ICF between the Continuum and Bowen+Hell data. A gaussian fit to the 
peak is over-plotted as the dotted-dashed line. The vertical lines mark the 
measured time-lag for each case: 11.8 ± 0.1s, 9.0 ± 0.7s and 2.2 ± 0.3s 
respectively 
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Figure 6. ICF for the selected 3min block of W3. The ICFs have been com- 
puted using the X-ray data and a B+Hell light curve obtained by subtracting 
from the green channel the Continuum data scaled by a factor cf = (i. e. 
no subtraction), 0.4 and 0.8 (optimal subtraction). The dotted line shows the 
3<t confidence level and the dashed vertical lines time-lags of 11.8 ± 0.1s, 
13.3 ± 0.1s and 16.0 ± 0.3s respectively. 

of continuum to be subtracted in order to amplify efficiently the 
reprocessed signal from the companion. Table [2] presents the ICF 
parameters obtained after subtracting a fraction cf (with cf in 
the range 0-1) of the red continuum from the B+Hell light curve. 
The table lists both the delay and the correlation level (obtained 
through a Gaussian fit) as a function of cf. It is clear that the 
higher is cf (i. e. more continuum is subtracted) the lower the 
correlation level of the peak. On the other hand, the ICF peak 
smoothly shifts to longer delays for larger cf values, as expected 
if more disc contribution is subtracted. We find that for cf values 
greater than ~ 0.8 — 0.85 the ICF functions become noisier and 
secondary peaks start to appear with correlation levels comparable 
to the main peak. Since for cf — 0.80 the correlation level of the 
peak is still significant with the 99.99% confidence level (i. e. 4cr) 
we decided to take this value as our best estimate of the continuum 
subtraction level for the ICF analysis. We call this 'optimally' 
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Table 2. Time lag and correlation level versus Continuum subtraction factor 



Continuum factor (cf) 


Correlation level 


Mean delay±lcr(s) 


0.00 


0.48 


11.8 ±0.1 


0.10 


0.45 


12.1 ±0.1 


0.20 


0.42 


12.4 ±0.1 


0.30 


0.39 


12.8 ±0.1 


0.40 


0.35 


13.3 ±0.1 


0.50 


0.31 


13.8 ±0.1 


0.60 


0.26 


14.5 ±0.2 


0.70 


0.22 


15.4 ±0.2 


0.80 


0.18 


16.0 ±0.3 


0.85 


0.15 


16.5 ±0.4 



subtracted lightcurve 'B+Hell CS', although it should be noted 
that this definition is somewhat arbitrary and could depend on 
our signal-to-noise. We also note that the la errors quoted in the 
table are formal and they are unrealistically small. Therefore we 
decided to choose a more conservative delay in the range 14.3- 
16.3s, corresponding to c/=0.6-0.8, when the la errors start to rise. 

In order to best estimate the actual amount of subtracted 
continuum it is necessary to take into account both the different 
throughputs of the filters and the quantum efficiency of the detec- 
tor, together with the spectral distribution of Sco X-l. For the case 
of a flat fx distribution between ~ 4600 A and 6000 A, we esti- 
mate that ~ 95% of the continuum level has been subtracted. On 
the other hand, i f we assume the spe ctral distribution reported by 
ISchachter et al.l <f 1989; hereafter S89) the actual subtraction level 
for cf — 0.8 would be ~ 60 — 70%. However, it is clear that the 
spectral distribution must depend on the X-ray state and the orbital 
phase, and therefore simultaneous optical spectroscopy would be 
needed to constrain the continuum level at the Bowen blend spec- 
tral region at the time of the observations. Note that only the nar- 
row components of the Bowen Blend arise from the companion, 
and hence the disc contribution can not be totally removed by sub- 
tracting continuum light from the B+Hell emission. Since the de- 
lay obtained for the B+Hell CS is ~ 3 — 4s higher than for the 
non subtracted data, this analysis strongly suggests that the com- 
panion signal dominates over the disc emission when most of the 
continuum contribution is removed from the B+Hell data. This is 
illustrated in fig.|6]where we present the ICFs for c/=0, 0.4 and 0.8. 
The plot clearly shows the monotonic shift of the peaks to longer 
delays when cf increases. Note that we obtain similar results by 
fitting parabolic functions to the ICF peaks instead of Gaussians. 

3.1.3 Fitting transfer functions 

In a second step we decided to fit the ULTRACAM light-curves 
with the results of convolving the RXTE data with a set of synthetic 
transfer functions. For simplicity we have used Gaussian TFs with 
two free parameters which represent the mean delay (to) and the 
standard deviation of the delay (At): 

TF(t) = Ae^s^i (2) 

where A is a normalization constant. This TF, which was first used 
in H98, is different to that computed by MD05 for the specific case 
of reprocessing in emission lines. However, we think that this is 
more appropriate since the data quality is not good enough to jus- 
tify the use of such an accurate model which considers separately 
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Figure 7. Gaussian Transfer Function fit for the SDSS v! , Continuum, 
B+Hell and B+Hell CS data. We show the la and 2.6<r contours for all 
the cases. Both TF parameters, to and Ato, are highly constrained by the 
fit. 
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Figure 8. X-ray light-curve (upper panel) and its convolution with 
TF(x^ in ) superimposed on the B+Hell (middle panel) and B+Hell CS 
(lower panel) data. In both cases the convolution between the X-ray data 
and TF{x 2 min ) for r = is over-plotted as a dashed line. The light- 
curves have been smoothed for clarity. 



the companion and disc contributions. We have applied this tech- 
nique to all data sets (i.e. B+Hell, B+Hell CS, SDSS it' band and 
Continuum) and find to values completely consistent with those 
obtained in the cross-correlation analysis. Fig.[7]shows the minima 
of the xt distribution for all the filters together with the la confi- 
dence level. A summary of the results is presented in Table[5] 



Table 3. TF and ICF analysis for W3 





ICF±1<t(s) 


r ± 1<t(s) 


At ± lc(s) 


B+Hell CS 


14.3-16.3 


13.5 ±3.0 


8.5 ±3.5 


B+Hell 


11.8 ±0.1 


10.75 ± 1.25 


6.25 ± 1.5 


SDSS u' 


9.5 ±0.7 


9.0 ± 1.5 


5.5 ±2.0 


Continuum 


9.0 ±0.7 


8.5 ± 1.0 


5.5 ± 1.5 



This new independent analysis confirms the time-lags ob- 
tained with the ICF method for the three channels. We have also 
applied this technique to the continuum subtracted data as we 
did for the ICF study. Again, we find that the mean delays shift 
to higher values when the continuum is subtracted. However, for 
cf > 0.5 — 0.6 the la contours become considerably large as well. 
Although the maximum delay is obtained for cf — 0.8 we have se- 
lected cf — 0.7 as our preferred subtraction level for the TF anal- 
ysis since the lower limit of the la contour hits its maximum value 
and hence is more constraining. This results in to = 13.5 ± 3.0. 
On the other hand, the Ato values are in the range 5-7 seconds for 
the three optical bands, whereas for the subtracted light-curve we 
get Ato = 8.5 ± 3.5s. Since the TF fitting is based on a Gaussian 
smoothing of the X-ray light-curve, the method tends to remove the 
uncorrelated features by increasing Ato and then does not attempt 
to fit the high frequency variability. However the fits successfully 
reproduce the main behaviour of the optical light-curves. We show 
in Fig. \E\ the fits for the optical (B+Hell and B+Hell CS) data for 
TF(Xmin)- We also show (dashed lines) that to = is clearly 
inconsistent with the optical observations. 

3.2 Correlated variability during W5 

We also find evidence for correlated variability during the last 
RXTE window on the night of 18 May (W5). W5 is centered at 
4>orb = 0.66 and during the observation the X-ray flux remained at 
6000 counts s _1 (i. e. half of the count rate in W3). In figure[9]we 
show a 15 min block of data during W5 where the correlated vari- 
ability is clearly detected. Following Sect. 3.1, we have analysed 
the data using both ICF and TF. 

3.2.1 ICF analysis 

The light curve obtained during W5 (see fig. |9) is essentially flat 
and only high frequency (> 0.002 Hz) variability is present. For 
the ICF analysis we have considered several data intervals within 
W5 and finally selected a ~ 6 min block which corresponds to the 
first part of the data presented in figure [9] (from day no. 4.202 to 
4.206) and the correlation is best seen. Note that the same delays 
are obtained by using the entire W5 data set, although for this 6 min 
segment the correlation peaks are more symmetric and hence the 
delays are better extracted from a gaussian fit. The ICF obtained for 
the three ULTRACAM bands (i.e. SDSS u , B+Hell and Contin- 
uum) clearly show highly significant peaks centered at 5.5 ± 0.2s, 
7.6. ± 0.4s and 9.4 ± 0.4s respectively as obtained through Gaus- 
sian fits. The ICF function for the B+Hell case is shown in figure 
[TUl The shapes of the ICFs are very similar for the three bands and 
have a peak correlation level of ~ 0.25. This value is much higher 
than the 3a confidence level obtained for this data interval, which 
is at 0.04. 

In order to obtain light-curves of the emission lines we have sub- 
tracted the red (Continuum) from the green(B+HeII) channels in 
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Figure 9. 1000s of X-ray (above), B+HeII (middle) and Continuum (below) 
data during W5. 
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Figure 10. ICF obtained during W5 for the X-ray and B+HeII data. The 
dotted line shows the 3<r confidence level and the dashed one the resulting 
time-lag of 7.6 ± 0.4s. 



Table 4. TF and ICF analysis for W5 





ICF±l<r(s) 


T ± 1ct(s) 


At ± lcr(s) 


B+HeII 


7.6 ±0.4 


8.5 ±1.5 


7.0 ±2.0 


SDSS u' 


5.5 ±0.2 


7.5 ±3.5 


18.5 ±4.0 


Continuum 


9.4 ±0.4 


9.25 ± 1.75 


10.5 ± 2.5 





Bowen + Hell / 
I 

Continuum 



-T:-*.- 



Figure 11. Gaussian Transfer Function fit for both B+HeII and Continuum 
data from W5. The dashed lines mark the to found for each case. The con- 
tours give the la confidence level for tq and Atq. 



in the X-ray data which are not present in the optical light-curves 
(see section 5.). We also obtain Atq = 7.0 ± 7.0s for B+HeII and 
this is consistent with the results for W3 within the errors. Figure 
[12] shows the fits of the convolved X-ray light-curves overplotted 
on the optical data. Although there are small features which are not 
correlated, the fit successfully reproduces the main behaviour of the 
optical data. 



the same way as we did for W3. We find that the peak's correla- 
tion level also decreases when cf increases but no shift to higher 
delays is observed in this case. In contrast to the results for W3, 
the ICF peaks tend to slightly shorter delays when the continuum 
is subtracted. This is consistent with the time-lags presented above 
for W5. 

3.2.2 TF analysis 

In the same way as we did for W3 we have computed synthetic 
optical light-curves by convolving the X-ray data with a Gaussian 
TF (see eq. [2}. Figure QT] presents the \ 2 ma P obtained from the 
entire W5 window(see fig.[9]l. We obtained to = 8.5 ± 1.5s and 
to = 9.25 ± 1.75s for the B+HeII and Continuum data respec- 
tively. These results are consistent with those obtained with the ICF 
method. In table [4] we summarize the results obtained with both 
techniques. 

Although the fit is poor we obtain t = 7.25 ± 3.5s for the u' 
band which is consistent with the ICF analysis. On the other hand, 
we get At = 18.5 ± 4.0s and At = 10.5 ± 2.5 for the u and 
Continuum bands respectively. These values are larger than those 
obtained for W3 and this is probably due to the presence of features 



4 COLOUR-COLOUR DIAGRAM ANALYSIS 

The broad band X-ray spectral variations of Z-sources 
are usually studied by building colour-colour diagrams 




200 400 600 800 

TIME(s) 



Figure 12. X-ray emission (top panel) convolved with TF(x^ n in) 
plotted on top of B+HeII (middle panel) and Continuum(lower panel) data. 
The light-curves have been smoothed for clarity. 
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Figure 13. Upper panel: colour-colour diagram obtained for the whole ob- 
serving campaign. During the first and third night the source was in the 
NB whereas it moved to the FB during the second night, where correlated 
variability is found. In the lower panel we separate the data of W3 and W5 
finding that Sco X- 1 was in a different X-ray state in the two epochs where 
we detect the echoes. 



( Hasin ger & van der Ml [1989). We have used the PCA en- 
ergy bands 6(2-4 keV), c(4-9 keV) and d(9-20 keV) in order to 
compute the Hard(d/c) and Soft(c/fc) X-ray colours. The colour- 
colour diagram obtained from the data of the whole campaign 
is presented in the upper panel of figure 1131 Sco X-l was at the 
bottom of the Normal Branch (NB) during most of the 17 May 
data but moved towards the Flaring Branch (FB) in the last RXTE 
visit, showing a 50 percent increase in flux. On 18 May Sco X-l 
was in the FB and exhibited large amplitude variability, with large 
flares similar to those seen during W3. On 19 May the system 
stayed on the NB again. As we showed in the light curve analysis, 
we have only detected correlated variability during the second 
night when the system moved across the FB. In the lower panel of 
fig. Q~3]we present the colour-colour diagram for the night of 18 
May and separate the W3 and W5 data using different symbols. 
We have also marked the intervals within these windows where 
the ICF analysis has been performed. It is clear that the ICF peak 
obtained during W3 is associated with the hardest X-ray emission, 
whereas during W5 (where lower delays are detected for the 
B+Hell emission) the system was close to the transition between 
the FB and NB. This study suggests that reprocessing associated 
with the companion star is most easily detected when the X-ray 
flux is at maximum and the spectra are harder. Note that although 



the spectra is harder in W3, the flux corresponding to both the soft 
and hard X-ray emission components is at maximum during this 
window. 



5 DISCUSSION 

We have found clear evidence for optical echoes to the X-ray 
emission when Sco X-l is moving across the FB. We have detected 
two different delay time scales by using both the low and high 
frequency variations present in the both X-ray and optical data. 
The low frequency analysis results in large delays of ~ 40 — 70s 
which are clearly not consistent with the RTS expected for Sco 
X-l if we only assume light-travel times. Such large delays have 
been previously seen by other aut hors. In particular a ~ 30 s delay 
was reported by 180 in Sco X-l. ICominskv et 1111(1987)1 studied 
the local diffusion time-scale during X-ray reprocessing. They 
found that 50% of the X-ray photons are reemitted within the first 
~ 0.2s, although a long tail of delays up to 10s is found as well. 
Since the light travel times within Sco X-l should be lower than 
15- 18s, it is clear that the time-lags obtained for the low frequency 
variability are not compatible with X-ray reprocessing on any 
possible site within the binary. Although the nature of these long 
time-lags is unclear we note that they are marginally consistent 
with the time scale for re-adjustment to thermal equilibrium in the 
accretion disc, which is exp ected to be about a few minutes (e. g. 
iFrank, King, & Raise! 1993) . 



On the other hand, the higher frequency variability component 
yields shorter delays (10 — lis) which are indeed consistent with 
X-ray reprocessing within the binary. In particular, during W3 we 
have found that the SDSS u band, the B+Hell emission and the 
Continuum data lag the X-ray emission with light travel times of 
~ 9s, ~ lis and ~ 9s respectively. In order to measure these 
delays we have used both ICF and TF methods finding excellent 
agreement between the two techniques. The case of the B+Hell 
emission is especially interesting since SC02 demonstrated that it 
is partially associated with emission arising from the companion 
star in Sco X-l. This emission is produced by fluorescence 
resonance through cascade recombination of Hell Lya photons 
in a low density medium, and hence we expect almost instanta- 
neous reprocessing times. Fig. [14] shows the range of expected 
time-delays for both the companion and the accretion disc as a 
function of the mass ratio (q — Mi/ Mi). In this simulation we 
have considered the TF presented in MD05 for the standard orbital 
parameters (i = 50°, Mi = 1.4M ) of Sco X-l and (f> orb = 0.52. 
The time-lags of ~ 10 — lis obtained for the B+Hell (unsub- 
tracted) data during W3 lies at the lower edge of the companion 
region, and shows that part of this emission must arise from 
reprocessing on the donor star. Hence, this time-delay represents 
the first direct detection of delayed echoes from the companion star. 

In an attempt to isolate the emission from the B+Hell lines we 
subtracted the Continuum data from the green channel and we have 
applied both the ICF and TF techniques to the obtained light-curve. 
By using the TF technique we obtain a time-lag of 13.5 ± 3.0s 
which is perfectly accommodated in the companion region (fig. 
I14t . On the other hand, the ICF method yields a time-lag between 
14.3 — 16.3s (see table |2). According to the TFs computed in 
MD05, this time-lag is larger than expected for the maximum re- 
sponse of the companion for each value of q (white dashed line). 
Here, we have assumed a disc flaring angle of a — 6° but the dis- 



Echoes from the companion 9 



15 



5 









Inclination: b0° _ 
■ ■ ■ ■: :■ . ■ 












BOWEN-Hcll |Bj 




CONTINUUM H 




DISK 



0.0 0.2 0.4 0.6 0.8 

MASS RATIO (q=M 2 /M 1 ) 



Figure 14. Expected range of delays for both the companion and the accre- 
tion disc considering i = 50°, Mi = 1.4Mq and the orbital parameters 
of Sco X-l. The black region shows the 1<t delays obtained for the B+Hell, 
and Continuum data. The shadowed region represents the range of delays 
obtained by subtracting the Continuum from the B+Hell light-curve. The 
white dashed line is the highest probability delay for the companion ac- 
cording to MD05 and using a = 6°. 



crepancy is still present if we use larger a values. This suggests 
that some of the values used for the orbital parameters of Sco X-l 
might not be correct. Furthermore, if we con sider that the disc ex- 
tends up to the truncation radius computed in |Paczvnskill(1977)l 
find that the time-delay of ~ 9s obtained for the Continuum is also 
slightly larger than expected for the outer disc. But this difference 
could be reconciled if the Continuum light were a combination of 
both reprocessed emission from the companion and the accretion 
disc. This is supported by the orbital modulation seen in the opti- 
cal light-curve, which is associated with heating of the companion 
(e. g. SC02). However this is in contradiction with the ~ 9s lag 
also measured for this filter during W5. Since W5 is centered at 
<t>orb = 0.66 the delay associated with reprocessing on the com- 
panion must be lower than for W3 (4> OT b = 0.52), and hence the 
measured time-lag should also be lower. We propose two possible 
scenarios in order to better accommodate our measured ICF time- 
delays with the model computed for Fig.l 141 
(i) The inclination angle of the system is ~ 60°, significantly 
higher than the i = 44 ± 6° reported by F01. This latter result 
was obtained by assuming that the jets are perpendicular to the or- 
bital plane of the binary. This assumption is difficult to test b ut there 
is evi dence for misalignments in some X-ray binaries (see iFendej 
120031 and references therein). 

(ii) The neutron star (NS) has a mass of ~ 1.8M©. Note that there 
are evidents of heavier (Mns > I.6M0) n eutron stars in some 
LMXBs (e. g. Cy g X-2 ;CCK98, V395 Car; IShahbaz et alj|2004 
llonker et aljBoO^ or XI 822-37 1:MCM051 

Obviously, a combination of the two scenarios can also solve 
the discrepancy i.e. i ~ 55° and NS mass ~ 1.6M@. We stress 
that this interpretation relies on the result derived through the 
continuum subtraction analysis for the ICF method. Alternatively, 
we obtain a more conservative result by using the TF method 
which results in time lags consistent with the standard parameters 
of Sco X-l. 

The ~ 9s delay obtained for the SDSS u band during W3 is 
much larger than measured during W5 (~ 6s). This is similar to the 
behaviour of the B+Hell emission where a ~ lis lag was measured 



during W3 and ~ 8s during W5, but different to the Continuum fil- 
ter, where a ~ 9s delay is observed both in W3 and W5. SC02 
have detected narrow components within the Bowen blend which 
arise from the companion star and also have reported the presence 
of many other high ionization lines at A < 4700 in the spectrum 
of Sco X-l. They showed that at least one of these emission lines 
also arises from the companion star. According to the UV spectra 
of Sco X-l presented in S89 there are many high excitation lines 
within the spectral range covered by the SDSS v! band. In partic- 
ular, this range covers three Bowen OIII lines at AA3030 — 3750, 
which probably also arise from the companion. The contribution of 
these lines to the u band light could explain the larger time-delays 
obtained for this band during W3. 

IWillis et all 1(1980; hereafter W80)l pointed out that the intensity 
of the high excitation emission lines present in the UV spectrum 
(AA1100 - 1800) of Sco X-l increases by a factor 1.5-3.2 when 
the source moves towards a high UV state. These authors associ- 
ated this high UV state with a high energy X-ray stat e by using 
simultaneous UV/X-ray data (see also lWhite et aljl976l) . Based on 
this study we propose that while Sco X-l stays at the bottom of 
the FB (W5) the intensity of the high excitation lines is not large 
enough to dominate over the continuum signal which arises from 
the accretion disc. And we note that the distribution of time-lags is 
correlated with wavelength for W5. This is consistent with a sce- 
nario where optical emission arises through reprocessing in an ac- 
cretion disc with a radial temperature profile. Moreover, when we 
subtract the Continuum from the B+Hell data in W5 the delays do 
not shift to higher values. This shows that the B+Hell emission is 
dominated by the same signal as the Continuum during W5. 
On the other hand, W3 is centered at orbital phase ~ 0.5 when 
the inner hemisphere of the companion presents its largest visibil- 
ity with respect to us. Our colour-colour diagram shows that at this 
time the source was at the top of the FB, which resulted in higher 
X-ray emission (factor ~ 2 with respect to W5) and a harder X-ray 
spectrum. According to W80, during this high energy state the in- 
tensity of the Bowen emission lines could have increased enough to 
dominate over the continuum signal and yielding the echo detection 
obtained for the W3 B+Hell data. On the other hand, we note that 
during the first part of W3 the system was also at the top of the FB 
but no time-lag consistent with the companion was detected. This 
shows that there are still some gaps in our knowledge of when and 
how it is possible to detect the reprocessing signal from the com- 
panion. 

For the case of the continuum subtracted data, where the errors 
are larger, we find that the accuracy provided by the ICF method 
is much higher than obtained through the TF technique. Neverthe- 
less, the latter has proved very useful to determine time-lags and 
their corresponding un certainties for the three bands considered. 
llvlcGowanetail [(2003) applied this technique to the sets of Sco 
X-l data used by 180 and P81. They measured r = 8.0 ± 0.8s 
and A ro = 8.6 ± 1.3s for the Johnson B band data. This filter 
covers our B+Hell narrow band and these results are consistent 
with ours during W5. The TF method also provides information 
about the time-delay distribution, which is encoded in the parame- 
ter Atq. During W3 we have measured Ato ~ 5 — 6s in the three 
bands whereas it increases after subtracting the Continuum from 
the B+Hell emission. These results are consistent with the time-lag 
distribution expected for reprocessing within the binary. Although 
the increase in Ato detected for the B+Hell CS data is consistent 
with the higher time-lag measured for this light-curve, we note that 
Ato is also sensitive to the presence of uncorrelated features in the 
X-ray data. Since the TF fitting is based on a Gaussian smooth- 
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ing of the X-ray light-curve, the method tends to remove the non- 
correlated features by increasing Atq. This clearly happens in the 
TF analysis of the SDSS u band during W5 where uncorrelated 
features are present. 



6 CONCLUSION 

We, for the first time, have isolated emission line contributions 
from the mass donor star that respond to X-ray variations with a 
delay that is consistent with light-travel time bewteen the X-ray 
source near the neutron star and the irradiated face of the donor 
star that directly contributes to the emission. The correlated signal 
from the companion is found at an orbital phase close to the 
superior conjunction of the donor, when a large flare occurs and 
the system stays at the top of the FB. In particular, the B+HeII 
emission is dominated by the reprocessed emission from the 
companion which increases the time-lag by 25% in comparison to 
that obtained in a featureless continuum at \ e ff = 6000. After 
subtracting the continuum emission from the B+HeII light-curve 
we amplify the signal from the donor and measure time lags in the 
range 10.5-16.5s by using both the TF and ICF techniques. We 
also note that for the latter case we have detected 4a significant 
ICF peaks centered at 14- 16s which suggests that either the orbital 
inclination or the NS mass could be higher than the standard values 
assumed. We stress here that a second detection of a B+HeII 
echo from the companion at a different orbital phase would 
directly constrain the inclination angle of the system. A more 
detailed study using Hell Lya and Bowen line light-curves could 
probably increase the signal from the companion, thereby allowing 
a more accurate determination of the system parameters in Sco X- 1 . 
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